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I.  INTRODUCTION 

Finite  array  analysis  is  useful  both  for  determinating  the 
characteristics  of  finite  size  arrays  and  also  for  determining  edge 
effects  for  elements  close  to  or  at  the  edge  of  very  large  arrays. 

In  this  report  the  solution  to  the  finite  array  problem  is  expressed 
in  terms  of  an  aperture  half-space  admittance  matrix  and  coupling 
coefficients  (scattering  parameters)  between  elements  of  the  array. 

The  problem  of  calculating  the  mutual  coupling  between  aper- 
tures in  a perfectly  conducting  ground  plane  previously  has  been 
studied  by  other  authors  [1-7].  Borgiottl  [1]  obtained  an  expression 
for  the  mutual  admittance  between  two  Identical  radiating  apertures 
in  the  form  of  a Fourier  transform  of  a function  strictly  related 
to  the  power  radiation  pattern  of  the  element. 

In  the  first  of  two  papers  [2-3],  Mailloux  found  the  near  field 
coupling  between  two  collinear  open-ended  waveguide  slots  by  formu- 
lating the  problem  as  a set  of  simultaneous  Integral  equations  and 
solving  the  resulting  equations  approximately  by  expanding  the  aper- 
ture field  in  a Fourier  series.  In  the  second  paper,  Mailloux  found 
the  near  field  coupling  between  two  closely  spaced  open-ended  square 
waveguide  slots  by  a first-order  analysis  which  is  based  on  the 
method  of  moments  using  a single-mode  approximation  to  the  aperture 
field.  He  also  gave  an  improved  first-order  analysis  which  uses  a 
higher  order  mode  solution. 
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A paper  by  Cha  and  Hsiao  [4]  and  dissertation  by  Hidayet  [5] 
treated  the  finite  rectangular  waveguide  array  problem,  while  papers 
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by  SCeyskal  [6]  and  Bailey  and  BosCian  [7]  dealt  with  finite  circular 
waveguide  arrays.  Steyskal's  method  Is  similar  to  the  one  used  In  this 
report.  He  found  the  scattering  matrix  for  a finite  circular  waveguide 
array  by  using  two  dominant  orthogonal  modes  as  expansion  functions  In 
a method  of  moments  solution. 

Formulas  for  the  computation  of  the  aperture  half-space  admit- 
tance and  scattering  matrices  of  a finite  planar  rectangular  waveguide 
array  are  derived  In  this  report.  The  procedure  assumes  a cosine 
aperture  electric  field.  For  apertures  that  are  close  together  the 
normal  quadruple  half-space  admittance  integral  Is  analytically  reduced 
to  single  Integrals  which  are  then  evaluated  numerically.  For  apertures 
farther  apart,  the  quadruple  half-space  admittance  Integral  is  analytically 
reduced  to  double  Integrals  which  are  evaluated  numerically. 


I The  general  theory  of  solution  is  an  extension  of  that  presented 

in  a previous  paper  [8] . The  basic  procedure  is  an  application  of  the 

i 

I method  of  moments  to  an  integral  equation  formulation  of  the  problem. 

The  unknowns  to  be  determined  are  the  equivalent  magnetic  current  coef- 

I 

I flclents.  The  magnetic  current  is  equal  to  the  tangential  electric 

field  in  the  aperture  regions,  rotated  90°.  The  Topelitz  property  of 
the  admittance  matrices  for  the  linear  and  rectangular  lattice  arrays 
is  utilized  when  solving  for  the  unknown  magnetic  current  coefficients. 
The  computer  program  which  uses  the  solution  derived  In  this  report 
will  be  given  in  a subsequent  report  [15] . 
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II.  STATEMENT  OF  THE  PROBLEM 

Figure  1 shows  the  problem  to  be  considered  and  defines  the 
coordinates  and  parameters  to  be  used.  The  infinitely  conducting 
plate  covers  the  entire  z = 0 plane  except  for  the  apertures  which 
are  rectangular  in  shape  with  side  lengths  a'  and  b'  in  the  x and  y 
directions,  respectively.  The  feeding  waveguides  have  inside  side 
lengths  of  a and  b in  the  x and  y directions,  respectively.  Note 
that  all  feeding  waveguides  have  the  same  dimensions  (a,b)  and  all 
the  apertures  have  the  same  dimensions  (a',b').  Also  z is  less  than 
zero  in  the  waveguide  region  and  z is  greater  than  zero  in  the  half- 
space region.  The  excitations  of  the  waveguides  are  sources  which 
produce  the  dominant  mode  wave  traveling  toward  the  apertures.  The 
waveguides  are  assumed  terminated  in  matched  loads  when  viewed  from 
the  aperture. 

We  first  consider  the  problem  of  a single  waveguide-fed  aper- 
ture radiating  into  a half-space  region  (see  Fig.  2a).  The  equivalence 
principle  [9,  Sec.  3-5]  is  used  to  divide  this  problem  into  two  separate 
regions  as  follows  (see  Fig.  2b).  The  aperture  is  covered  by  an  elec- 
tric conductor.  The  fields  in  the  waveguide  region  are  produced  by 
the  impressed  sources  the  equivalent  magnetic  current  M 

M - n X ^ (1) 

over  the  aperture  region  with  the  aperture  covered  by  an  electric 
conductor.  The  fields  in  the  half-space  region  are  produced  by  the 
equivalent  magnetic  current,  -M,  with  the  aperture  covered  by  an 
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FEEDING  ^APERTURE 
WAVEGUIDE 


PLANE 


1.  Waveguide-fed  rectangular  apertures  in  a perfectly 
conducting  plane. 
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(a)  ORIGINAL  PROBLEM 
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CONDUCTORS 


HALF-SPACE 
— REGION 


CONDUCTOR 


(b)  EQUIVALENT  PROBLEM 


Fig.  2.  A single  waveguide-fed  aperture  radiating  into  half-space 
bounded  by  an  electric  conductor. 
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electric  conductor.  The  condition  that  the  equivalent  magnetic  current 
In  the  waveguide  region  is  +M  and  in  the  half-space  region  -M  ensures 
that  the  tangential  component  of  electric  field  is  continuous  across 
the  aperture. 

Another  necessary  boundary  condition  is  the  continuity  of  the 
tangential  component  of  magnetic  field  across  the  aperture.  The  tan- 
gential magnetic  field  over  the  aperture  on  the  waveguide  side, 


is  equal  to 


^t  -'t  — 


where 


" is  the  tangential  magnetic  field  due  to  impressed 
sources 

VC 

(M)  is  the  tangential  magnetic  field  due  to  the 
equivalent  magnetic  source  M. 

On  the  half-space  side  of  the  aperture  we  have 

hJ"  • hJ"(-M)  = - hJ®(M)  . (3 

The  last  equality  in  (3)  is  a consequence  of  the  linearity  of  the 
operator.  Note  that  H^"***,  and  are  all  computed 

with  an  electric  conductor  covering  the  aperture.  The  true  solution 
is  obtained  when  of  (2)  equals  of  (3)  or 


+ H,  (M)  - - St  . 

This  the  basic  operator  equation  for  determining  the  equivalent 
magnetic  current  M. 
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Although  (A)  was  derived  for  a single  waveguide  in  a ground 
plane,  it  applies  equally  well  for  an  N-element  waveguide  array.  In 
this  case  the  waveguide  region  is  considered  to  include  all  of  the 
waveguides  and  the  half-space  region  is  as  before.  For  the  multiple 
aperture  case,  (4)  becomes 


^‘^i  ~ j=r‘^i  ^ '“'^i 


i = 1,2,.. .,N. 

In  (5),  is  the  equivalent  magnetic  current  for  the  jth  aperture. 

The  subscript  i denotes  magnetic  field  evaluation  in  the  1th  aperture. 


= V,  M,  ( 

where  is  a complex  constant  to  be  determined  and  is  an  expan- 
sion function  to  be  specified.  Substituting  (6)  into  (5),  we  obtain 


V.  + I V (M^)  = - (?: 

1 i j ^ 

i = 1,2, ...,N. 

Next,  define  the  symmetric  product  <A,B>  of  two  vectors  A and  B by 

<A,B>  = |j  A • B ds  (8] 

where  the  Integral  is  over  all  aperture  regions.  Also  define  a set 
of  testing  functions  {W^»  1“1,2,...,n}  which  may  or  may  not  be  equal 
to  the  expansion  functions.  Then,  taking  the  symmetric  product  of  (7) 


with  the  testing  funclon  W , , we  obtain 

—a 
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N 

V.  <W  )>  + I V )>  » - <W.,H^P>  < 

^ ^ ''i  ^ j=l  J 1 «^i  J 1 

i = 1,2, ...,N. 

Solution  of  this  set  of  linear  equations  determines  the  coefficients 

V and,  therefore,  the  equivalent  magnetic  current  M'^ . 
j 

Equation  (9)  can  be  rewritten  in  matrix  notation  as  follows: 
Define  an  admittance  matrix  for  the  waveguide  regions  as 


<-W^, 

where  is  the  Kronecker-delta  function 


1 i =■•  j 


0 i J 


and  for  the  half-space  region  as 


[y;;]  - [<-w^, 

The  minus  signs  are  placed  in  (10)  and  (12)  on  the  basis  of  power 
considerations.  Define  a source  vector 


i tj  '1 


and  a coefficient  vector 


V ■ fV  1 
^ i^Nxl 


The  resulting  matrix  equation  which  is  equivalent  to  (9)  is 


[Y«g  + y‘’®]V  - . 


I 

I 
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The  physical  Interpretation  of  (15)  is  that  of  two  generalized 

WO  ho 

admittance  networks,  [Y  and  [Y  ] , in  parallel  with  the  current 
source  By  inverting  (15),  we  obtain  the  resulting  voltage  vector 

V which  is  the  vector  of  coefficients  which  determines  M 


V = [y''®  + Y^®]"^ 


(16) 


The  expansion  (M.)  and  testing  (W.)  functions  are  defined  as 
-*■1 


where 


iJi  ° iii  ■ ■ 


Pj^(x.y)  = 


- a'/2<x<x^  + a'/2 


L 


- bV2 1 y 1 y^  + b'/2 


V 0 all  other  x,y 


(17) 


(18) 


The  equivalent  magnetic  current  over  the  ith  aperture  region  is 


(19) 


The  electric  field  in  the  1th  aperture  is 


£i  (x-x^)  . 


(20) 


We  choose  the  minus  sign  in  the  magnetic  current  expression  (19)  to 
avoid  an  explicit  minus  sign  in  our  expression  for  aperture  electric 
field  (20). 
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III.  ADMITTANCE  FORMULATION 


a)  Determination  of  - 


To  evaluate  the  aperture  admittance  (10)  In  the  waveguide  region, 
we  consider  a single  expansion  function  M^  on  the  z > 0 plane  In  the 
1th  waveguide  region.  The  tangential  field  produced  by  ^ can  be  ex- 
pressed In  modal  form  as  [9,  Sec.  8-1] 


k 


- - I AikV  H.  “ Sk 

k 

where  A^^^  are  modal  amplitudes,  modal  propagation  constants, 

are  modal  characteristic  admittances,  and  ^ are  normalized  modal 
vectors.  The  modal  vector  orthogonality  relationship  Is 


If 

guide 


0 1 j 


1 i - j 


where  the  Integration  Is  over  the  waveguide  cross  section.  At  z ~ 0, 
we  have 

Ma  " u “ I A.,  u X e.  . (: 

-1  -z  --t  -1  - “ Ik  -««  -k 

z"0  k 

Multiply  each  side  of  this  equation  scalarly  by  x e ^ and  Inte- 
grate over  the  waveguide  cross  section  obtaining 


fl  Kl  • H,  ■>  Sj  a*  - J Aj^  II 


(u  X t ) • (u  X p ) ds  . (24) 

^ *-z 


guide 


guide 


By  orthogonality  (22),  all  terms  of  the  summation  are  zero  except  the 
j-k  term.  Hence 


I 


r 


A **■ 
*ik 


If  £i  • 1 


u X e,  ds  . 

-Z  "“IC 


apert. 


We  have  replaced  the  integral  over  the  waveguide  cross  section  by  one 
over  the  aperture  since  exists  only  in  the  aperture  region.  Sub- 
stituting the  second  equation  of  (21)  evaluated  at  z * 0 into  (10) , 
we  obtain  the  diagonal  elements  of  the  waveguide  admittance  matrix. 


If  Si-!!: 


X e,  ds  . 
z 


apert. 


Since  M,  * W , , (26)  becomes 
-►i  —1 


The  set  e.  of  modes  for  the  rectangular  waveguide  is  split 
— k 

into  a set  e?^  of  TE  modes  given  by  [9,  Equations  (8-34),  (3-86), 
—k 

and  (3-89)  and  Section  4-3] 


V (mb)2  + (na)2  ® 


m . mnx  mry  , 

- u — sin  cos  -r^  ] 

^ a a D 


in  * Xi 


n • L 


+ n j*  0 


1 m ■ 0 

2 m - 1,2, . . . 


and  a set  e^^  of  TM  modes  given  by 


.1 
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TM  - / ab  r n>  “'Tx  . niry 

e . / * 2 * / 5 TT  [u  — cos  sin 

V („b)2  ^ ^^^^2  -*  a a b 


, n . mTix  nTTv, 

+ u T sin cos  -r^] 

b a b ' 


m ■ 1,2,3, ..  .L 
m 

n * 1,2,3,...,L 


Note  that  equations  (28)  and  (29)  are  valid  only  when  the  origin  is  at 
a corner  of  the  waveguide. 

Substituting  (17),  (28),  and  (29)  into  (25),  we  obtain 


fl  £!i  • S 


u X ds 
--z  ~ic 


apert. 


mil  . / ran  ^ , \ 


where 


k ■ m+n(L  +1) 


m»0,l,2, . . . ,L 
n 

n*0 


m+n  ^ 0 


■ If  a 

apert. 


TM 

M.  • u X P ds 
*il  --z  -“k 


k^b^/aa’^b' 


k - iirf(n-l)L 


’’m-1,2,3, . . . ,L 
n 

“ 1 n-1,2,3 L 

^ n 
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I(m,n) 


(b+b’)/2  (a+a')/2 

II  dy  C08  II  dx  cos  (x-a/2)  sin 

(b-b')/2  (a-a^/2 


2b* 


Tra*  2 - 

(H L.) 

'■  2 ,2^ 

a a* 


. nTTb* 

^ mw  mr  mTra*  ® " 2b 

sin  -x-  cos  cos  — r tt 

2 2 2a  mrb 

2b 


(33) 


TE  TM 

Note  that  the  introduction  of  e.  and  e,  has  split  the  coefficients 

■"K  "tc  2 

(25)  into  aJ®  and  A™.  If  (a/a*  » m)  in  (33),  (cos  ^^/(%  - -^))  is 

2 a a* 

to  be  replaced  by  its  limit  (-ira  /4).  If  n is  zero  in  (33),  (sin(0)/0) 
is  to  be  replaced  by  unity. 


The  characteristic  admittances  of  a rectangular  waveguide  with 

relative  dielectric  constant  and  relative  permeability  unity  are 

TV  TM 

classified  as  either  TE  admittances  or  TM  admittances  Y^  given  by 


[9,  Section  4-3], 


,.TE 


r 


- ^ 


W'r-  'T> 


k^< 


kv€  > k 
r i 


(34) 


i ■ m+n(L  +1) 
tn 


in^O • • * y L 
11*0 


m 


m+n  ^ 0 


and 


r 


,TM 


1.  ' 


k»€  < k. 


V<T)  - S 


k*€^  > k^ 


1 ■ m+(n-l)L 


r tn"l,2,3,...,L^ 

“ I n-1,2,3, . . . ,L  . 

^ n 


(35) 
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In  (34)  and  (35),  n Is  the  characteristic  impedance  of  free  space. 


k is  the  free  space  wave  number,  and  is  the  cut-off  wave  number 


given  by  (32) 


b)  Determination  of  Y 


hs 


ij 


Since  the  apertures  are  covered  by  conductors,  the  z * 0 plane  is 
a complete  conducting  plane  and  image  theory  applies.  The  magnetic 
current  expansion  functions  are  on  the  surface  of  the  z ^ 0 plane. 


Their  Images  are  equal  to  them  and  are  also  on  the  z « 0 plane.  The 


,hs 


result  is  that  [Y  ] is  the  admittance  matrix  obtained  using  expansion 


functions  radiating  into  free  space  everywhere.  Therefore  (12) 


can  be  written  as 


,.h8 


.is 


- 2<W^,  (Mj)> 


(36) 


where  H (M. ) is  the  magnetic  field  in  the  1th  aperture  region  produced 
••t.  -“j 

^ .fs 


by  radiating  into  free  space.  The  magnetic  field  (M  ) can  be 

'“j 

expressed  in  terms  of  an  electric  vector  potential  F..  and  a magnetic 

■^ij 


scalar  potential  (|)^j  as  [lOJ 


.fs 


S.;  - 2 fy 


(37) 


where 


,,  -3k|£-S.’l 

£«  ■ jj  Sj  '|r  - r’l 

apert . 


(38) 


-Jkk-r’l 

♦ij  ■ J)  '|r  -r-l 

apert. 


(39) 


V • 

r 

-jO) 


(40) 


sdtmttsasaasuittEi 


! • 
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where  r and  r'  are  respectively  the  vectors  to  the  field  and  source 
points,  (1)  Is  the  angular  frequency,  e Is  the  permittivity  of  free 
space,  p Is  the  permeability  of  free  space, and  k Is  the  free  space 
propagation  constant.  Substituting  (37)  Into  (36)  and  using  (8),  we 


obtain 


yhs  _ 2 

ij 


**•1 

J J 


(jw  + V (J)^j)ds  . 


apert. 

Because  of  the  Identity 


f ff  f 

0 = V • ((()., W.)ds  = W.  • V((),  .ds  + V • W ds 

— ij**l  J J "1  " Ij  J J ij 


apert. 


(41)  becomes 


apert . 


apert . 


where 


JJ  ^^13  • «i  + <^ij 

apert. 


V • 

’l  “ - jO) 


Next  define  by 


^ If  hi  • Si  • 

apert. 


Substituting  (17)  Into  (45),  we  obtain 


y^+b'/2  x^+a'/2  yj+bV2  Xj+a'/2 

^IJ  * TTa^  I dy  I dx  cos  | dy’  | dx'cos^(x'-Xj)G(x'-x,y'-5 


yi-b'/2  x^-a'/2 


yj-b’/2 


where 


G(x'  - X,  y’  - y) 


'-x)^  + 


j-jkV(x'-x)‘  + (y*-y) 

V(x'-x)^  + (y'-y)^ 


Consider 


(x'-x) 


y^+bV2  yj+b'/2 

I dy  I dy'  G(x'-x,  y'-y)  . 
y^-b'/2  yj-bV2 


Substituting  the  transformation  y'  » v+y  Into  (48),  we  obtain 


l^j(x'-x) 


y^+b'/2  yj+b'/2-y 

I dy  I dv  G(x'-x,  v)  . 


"if  ' "j  " 
dy  d\ 

y^-b’/2  yj-bV2-y 


Interchanging  the  order  of  Integration,  we  obtain  (see  Fig.  3) 


(x'-x) 


yj-yj  y>bV2 

dv  G(x'-x,  v)  dy 


yryr*’’ 


yj-b'/2-v 


yj+b'/2-v 

I dv  G(x’-x,  v)  I dy 


yj-yi 


y^-b'/2 


Consider 


x,+a’/2 


X +a' /2 


- I dx  cos  ^ (x-x^)  I dx'  cos  ^(x ' -x ^ ) (x ' -x) 


x^-a'/2 


x^-a'/2 


Substituting  the  transformation  x'  - u+x  Into  (51),  we  obtain 


x>a'/2 


X .+a' /2-x 


fa.  f a.  A 

- I dx  cos  ^ (x-x^)  I du  cos  ■^(u+x-Xj)I^j  (u) 


x^-a'/2 


Xj-a'/2-x 


(52) 
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Interchanging  the  order  of  Integration,  we  obtain 


lx 

Ij 


I du 

Xj-x^-a' 


Xi+a'/2 


dx  cos  (x-x  ) cos  (u+x-x,) 

did  I 


Xj-a'/2-u 


x.-x.+a'  x.-ha’/2-u 

j fi  „ i r 

+ J du  lL(u) 


ij 


x^-aV2 


dx  cos  (x-x  ) cos  (u+x-x  ) , 

did  J 


Substituting  (50)  and  (53)  into  (46),  we  obtain 


1 _ jCDE  -lx 

IJ  TT  a’b'  ij 


Next  define  by 


Yy  - 21^ 


j 


ds. 


apert. 

Substituting  (17)  into  (55),  we  obtain 


-JTT 


y^+b'/2  x,+aV2 


toya ' ^b ' 


y,+b'/2  x,+a'/2 


(53) 


(54) 


(55) 


r ' ii-  ' ■'j  fjj  ' 

dy  j dx  sin  ^(x-Xj^)  j^y'J  sin  ^(x'-Xj)G(x'-x,y'-y) , 


y^-bV2  x^-aV2 


,j-b'/2  Xj-aV2 


(56) 


Substituting  the  transformations  x'  - u+x  and  y'  ■ v+y  into  (56)  and 
interchanging  the  order  of  integration,  we  obtain 


a iLH 

'J  u.ua*\- 


(57) 


where 


x,+aV2 


I.^  “ du  (u)  dx  sin  (x-x  ) sin  (u+x-x.) 
ij  J ij  J a ' i a j 


Xj-x^-a'  Xj-a'/2-u 


x.-x.+a' 

i ri 


x,+a'/2-u 


•i  j C 

du  (u)  dx  sin  (x-x  ) sin  (u+x-x  ) . 

ij  did  J 


Xj-Xi 


x^-a'/2 


After  a tedious  but  straightforward  evaluation  of  two  integrations  for 

lx  2x 

each  of  the  quadruple  Integrals  and  , we  obtain 

Ij  Ij  ij 


fj  i TTU  TTU 

7^7  { J dv(K^+v)[  J du  G(u,v)[(K2+K^u)cos  ^ +(K^+K5u)8in 


’rj'"i 


yj-Vi-b'  xj-x^-a* 


x.-x.+a' 

jf  i 


J f 1 

+ j du  G(u,v) [ (K^-K^ujcos  ^ + (Kg-K^ujsin  ^]] 


Xj-Xi 


+ I dv(Kg-v)[  I du  G(u,v) [ (K2+K2u)co8  ^ + (K^+KjU)sin 


Yj-yi  xj-x^-a 


Xj-Xi+a' 


+ f du  G(u,v) [ (K^-K.u)co8  ^ + (Kg-K^u) sin  ^] ] } 
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where 


Ki  = y^-y.  + b' 


2 2 
*^2  * (1-  cos  (Xj-x^)]  (1+ y^)sin[^(Xj-x^)] 


k a'  ■’ 


2 _ , 2 
K4  - (1-  -f '2^  (x^-Xj+a')  sin[^  (x^-x^)]  - ^ (1+  (^3) 


S “ 72^^  (Xj-x^)] 

k a ■' 


H ■ yj  - n * »’ 


K,  • (1-  -2^)<«j-»i+*')  ^ <1+  -j-j)sln[i(Xj-Xj)l 


2 , 2 
Kg  “ (1 (Xj-x^)]  (1+ -2^)cos[^(Xj-x^)] . (67) 


Substituting  the  coordinate  transformations  u = p cos  6 and  v = p sin  6 
into  (59)  and  integrating  out  the  p variable,  we  obtain  (see  Fig.  4) 


ij  27rTia 


2 

^ {|  de[K^[K2L^(0)  + K3L3(0)  + K^L2(0)  + K5L^(0)] 


+ [K2L5(0)  + K3L2(0)  + K^L^(0)  + K3Lg(0)]] 


+ I d0tK3(K2L^(0)  - K3L3(0)  + KgL2(0)  - K3L^(0)] 


+ [K^L3(0)  - K3L2(0)  + KgLg(0)  - K5Lg(0)]] 
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I 

i 


I 


deEK^lK^L^Ce)  + K^L^O)  + K^L^O)  + KjL^CG)] 


(HI) 


- [K^L^ce) 


de[Kg[K^L^(9) 


®6 

(IV) 


+ K3L^(e)  + K^L^(e) 
- K3L3(0)  + KgL2(e) 


+ KgLgO)]] 

- K3L^(e)] 


where 


- [K^L3(0)  - K3L2(0)  + KgLg(0)  - K3Lg(0)]]} 

M (0)  - M (0)  M (0)  - M (0) 

-•i'®'  ■ ^ < m,(6)  - -M,(e)  ■ > 


(68) 


(69) 


r 
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L-,(0)  = sin  0 cos  0{-  ^ [ 


^ P2(0)M^(0)  - pJ(0)M2(0)  P2(e)M^(0)  - p^(0)M^(0) 


M5(0) 


M6(0) 


P2(0)M^(0)  - P^(0)M2(0)  P2(0)M3(0)  - p^(0)M^(0) 


M, (0)-M,(0)  M-(0)-M, (0) 

+ j I- i rr — 1> 


(M5(0))- 


(Mg(0))' 


(75) 


Lg(0) 


^ p^(0)M^(0)  - p2(0)M2(0)  P2(0)M3(0)  - P3(0)M^(0) 

sin  0 cos  0{-  2 t TO) + TO) 1 


P2(0)M^(0)  - P3(0)M2(0)  P2(0)M3(0)  - P3(0)M^(0)^ 

j [ _ _ ] 


(M3(0))' 


M (0)  - M (0)  M (0)  - M.(0) 

+ i — + ^ 1 } 


(Mg(0))' 


(76) 


(M5(0))- 


(Mg(0))- 
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The  variables  pj^(0)  and  pjCQ)  are  respectively  the  lower  and  upper  limits 
of  integration  for  the  p integration  in  (59)  after  the  transformations 
u * p cos  6 and  v = p sin  9 are  used.  They  are  dependent  on  the  0 vari- 
able and  the  Integration  subarea.  For  instance,  when  integrating  over 
region  IV  between  9^  and  9^  (see  Fig.  4),  Pj^(0)  = (Xj-x^)/cos  0 and 
P2(9)  = '’’)/sin  0. 


IV.  SINGULAR  POINTS  OF  Y 


The  integrand  of  (68)  has  removable  singularities  when  ka'  < it. 
When  either  |cos  0 - 2a' | < e (small  value)  or  |cos  0 + 2a' 1 < e,  we 
replace  equations  (69)-(76)  by  their  respective  limits.  For  an  example, 
consider  equation  (69), 

, M,(0)-M_(9)  M (0)  - M (9) 

2 ^ MJ9)  M,(9)  ^ 

5 o 

jp2(9)(^  cos  0 - k)  jPj^(0)(^  cos  0 - k) 

= _ i {6 


I 

(^  cos  0 - k) 

-jp2(0)(^  cos  0 + k) 

P P 


+ k)  cos  0 + k) 

- e 


“ _ • j • 

(-^  cos  0 + k) 

SL 

If  |-^  cos  0 - k|  < E,  then  by  retaining  the  first  three  terms 

d 

of  the  Taylor  series  for  the  exponentials  in  the  first  half  of  (69), 
we  obtain 

Lj(0)  - - ^ {1  + jP2(0)(^  cos  0 - k)  - ■|•tP2^6)(aT  cos  0 - k)]^ 

- 1 - jPj^(0)(-^  cos  0 - k)  + 9 " 

cos  0 - k 
a 

M (0)  - M,(0) 

_ .J 2 } . 

M^(9) 


F 
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If  we  take  the  limit  as  — cos  0 approaches  k for  the  first  half  of 

3. 

(83),  we  obtain 

L,(6)  = j (P  (0)  - p (0))  .^3(0)  ~ 

^ M,(0)  ’’ 

o 

If  cos  0 + k|  < e,  then  by  retaining  the  first  three  terms 

3 

of  the  Taylor  series  for  tho  exponentials  in  the  second  half  of  (69) , 
we  obtain 

M, (0)  - M,(0) 


(84) 


L^(6)  " - 2 M (0) 


•^^)  + I (P2(e)  - P^(0)). 


(85) 


obtain 


If  we  apply  the  preceding  procedure  to  equations  (70)-(76),  we 


1 ”3(6)  - M^(e) 
L^O)  = - 1 (p^(e)  - p^(0))  - i (^-^-  : 


(86) 


L2(6)  “ cos 


L^(0)  = cos 


1 , , . p (0)M  (0)-p  (0)M. (0)  . M-(0)-M,  (0) 

e {i(p^(9)-pj(9))  + i(-^ „ re?  ) + ¥:~J  » <"> 

b (Mg(0)) 

.0  ? P„(0)M  (0)  -p  (0)M  (0)  M_(0)-M,(0) 

e {-  i(p2(e)-p2(e))  -i(-^ ^)  ^(-1--^)}  (88) 


(Mg(0))' 


^(0)  = 1^(0)  tan  0 


L^(0) 


L2(e) 


L^(0)  tan  0 

1 o , , Pn(e)M,(0)  - pJ(0)M,(0) 

sln0cos0  {|(pj(6)-pj(0))  -■  (Qp- ^ 

6 


(89) 

(90) 


) 


P2(0)M3(0)  - p,(0)M,,(0)  M,(0)  - M,.(0) 


(Mg(0))‘ 


j (- 


(Mg(0))' 


-)} 


(91) 


M 


. . o 1 p;(e)M-(0)-p;(0)M,(0) 

Lg(0)  » sin  0 cos  e {-  ^(PjCe)  - pjo))  - — ) 


p„(9)M,(e)  - P-(e)M  (6)  M (0)  - M,  (6) 

+ j (-^ — ^ 1 — —-i 

(Mg(0))  (Mg(0)) 


for  cos  0 - k]  < G,  and 


= - I ^ -M3(0)  — ^ 2 


Lj(6)  = cos 


. P_(0)M,(0)-p,(0)M-(0)  M (0)-M,(0)  , , , 

8 f-  i<  1(9)  ■ — > ^ 2 > + 


p (0)M  (0)-p  (0)M 
L,(0)  =cose{-  i( 


L^(0)  = Lj(0)  tan  0 


,(0)  . M.(0)-M-(0)  . , „ 

- — ) — V-)  +t(p;(0)-p7(0))}  (95) 

^ / r\  \ \ ^ ^ ^ i- 


(M5(0)) 


L^(0)  = L^(0)  tan  0 

p^(0)M  (0)-p,(0)M  (0)  p (0)M  (0)-P  (0)M  (0) 
L,(0)  - sin  0 cos  0 {-  ^ ^ ) + 1 

M,(0)-M,(0)  , - - 

+ j (— i ) + T (P,(0)  - P^(0))} 

(H^wy 

, p5(0)M, (0)-p?(0)M  (0)  p (0)M  (0)-p  (0)M  (0) 
L_(0)  - sin  0 cos  0 {-  STTfi^ ^ ^ ^ 2 ^ 

M,(0)-M,(0)  . , , 

+ i -*■  i (P,(0)  - p:(0))) 

/ A\  \ O ^ X 


(M5(0))- 


for  1^  cos  0 + k]  < e. 


(99) 


27 


v^__jcattert'‘’"  matrix  - [S] 


To  determine  the  scattering  parameters  5^,^  for  a waveguide  fed 


aperture  antenna  array  radiating  into  a half-space  region,  a source  is  placed 
in  waveguide  j and  all  waveguides  are  terminated  in  matched  loads. 
relates  the  backward  wave  traveling  toward  the  load  in  waveguide  1 to  a 
forward  wave  traveling  toward  the  aperture  in  waveguide  j . Besides  the 
field  directly  coupled  from  waveguide  j to  i,  there  are  components  due 
to  the  radiated  fields  from  the  other  apertures  resulting  from  induced 
aperture  fields.  Therefore,  the  scattering  matrix  explicitly  identifies 
the  Incident  and  reflected  (or  coupled)  wave  in  every  waveguide  element. 

Let  a source  which  excites  only  the  dominant  mode  be  placed  in 
the  jth  equivalent  waveguide  region  and  represent  the  fields  transverse 
to  the  2-direction  as 


, "^O^  ^O^  WE  i 

E (z)  = (e  - e ) e + 

wtj  -o  -t  - 


V 


(e  - e ) + Vj  I A^^^e  ^ 


(100) 


IL  (z) 


Y (e  ° + 

o 


e ) u 

^'7 


e 


''j  ^ \ ' 


u 


(101) 


where  the  index  o denotes  the  dominant  mode,  ate  the  modal  propagation 

constants.  A.,  are  the  modal  constants  defined  by  (25),  V,  is  the  unknown 
J ^ J 

magnetic  current  coefficient  determined  by  (15),  and  the  are  modal 
characteristic  admittances  defined  by  (3A)  and  (35).  The  electric  field 
transverse  to  the  z-directlon  in  the  ith  equivalent  waveguide  region  (no 
Impressed  source)  is 


1 
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k 


V 


h.  ■ 


(102) 


The  scattering  parameter  is  the  ratio  of  the  backward  traveling 
wave  in  waveguide  i to  the  forward  traveling  wave  excited  by  a source 
in  waveguide  j . (It  is  assumed  that  only  the  dominant  mode  wave  propa- 
gates in  the  waveguide  region.) 

For  i = j,  we  obtain  from  (100) 


‘ ''j  ^ 


(103) 


For  1 / j , we  obtain  from  (100)  and  (102) 


(104) 


To  determine  V^,  let  us  consider  equations  (15)  and  (13) 


where 


[Y 


(15) 


(13) 


is  equal  to  the  first  term  on  the  right  side  of  (101)  when  2*0. 

"‘'j 

When  the  aperture  is  covered  by  a conductor,  the  waveguide  is  terminated 
by  a conducting  plane.  According  to  image  theory,  the  tangential  magnetic 
field  at  2 » 0 is  then  just  twice  the  incident  wave,  or 


uimp 

-J 


2 Y u X e 
o -o 


(105) 


J 
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This  Is  the  used  in  (13)  to  evaluate  the  excitation  vector 
Hence,  the  components  of  the  excitation  vector  are 


limp  = 2 Y ff  W . u X e ds 
J o J J - j — z ^ 


apert. 


2 Y A, 
o jo 


ifP  = 0 


if  i - j 


if  1 j . 


(106) 


(107) 


For  an  equally  spaced  linear  array,  the  matrix  [y''®  + y''®]  in 
(15)  is  symmetric  Toeplltz.  Equation  (15)  is  solved  by  a recursion 
relation  given  by  Zohar  [11]  which  is  based  on  an  algorithm  by  Trench  [12]. 

For  a uniformly  spaced  rectangular  array,  the  matrix  [y”®  + Y*'®]  in 
(15)  is  symmetric  block-Toeplitz.  It  is  symmetric  since  we  have  used  a 
Galerkin  solution  (choosing  expansion  functions  equal  to  the  testing 
functions) . The  block-Toeplitz  property  for  a uniformly  spaced  rectangular 
array  is  given  by  a matrix  of  the  form 


Y Y-  . . 

o 1 


(108) 


Y Y , 
n n-1 


where  Y^  is  that  submatrix  which  defines  the  self-admittance  of  an  ele- 
ment of  the  array  and  Y|^_j|,  i j is  that  submatrix  which  defines  the 
mutual  admittance  between  elements  1 and  J of  the  array.  An  element 
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(submatrlx)  for  the  rectangular  array  is  either  the  shortest  row  or 
column  of  apertures  in  the  x-y  plane  (see  Fig.  5).  A recursion  rela- 
tion which  uses  this  block-Toeplitz  symmetry  property  [13]  is  used  to 
solve  (15) . 

For  uniformly  spaced  Isosceles  triangular  arrays,  equation  (15) 
is  solved  by  Gaussian  elimination  and  LU  decomposition  [14] . 

It  should  be  mentioned  that  all  three  methods  used  in  solving 
the  linear,  rectangular,  and  triangular  arrays  solve  (15)  directly 

W2  llS 

without  computing  the  inverse  matrix  to  [Y  ° + Y ] . 

VI.  SAMPLE  COMPUTATIONS 

A computer  program  using  the  formulas  derived  in  the  preceding 

sections  has  been  written.  It  is  described  and  listed  in  a subsequent 

report  [15].  In  this  section  we  give  some  examples  of  the  computations 

that  can  be  made  using  the  general  program. 

Figure  6 shows  the  half-space  admittance  between  two  narrow  slots 
hs 

in  echelon  [^2.2^  separation  distance  is  varied  while  the  angular 

position  is  held  fixed.  This  example  tested  the  limiting  expressions 
derived  in  Section  IV  since  for  0 ••  0,  (-^  cos  0 - k)  ■*  0.  Our  compu- 

d 

rations  are  compared  to  those  of  Borgiottl  [1]. 

Figure  7 shows  the  scattering  parameter  (amplitude  and  phase) 
for  two  slots  as  the  angular  position  of  one  is  varied  with  respect  to 
the  other  while  the  separation  distance  is  held  fixed.  Our  computations 
are  compared  to  those  of  Mallloux  [3]. 


Fig.  5.  Choice  of  elements  for  a block-Toeplitz 
admittance  matrix  for  a uniformly  spaced 
2x3  rectangular  aperture  array. 


MHOS  X 10-* 


32 


SA 


Fig.  6.  The  half-space  admittance  between  two  narrow  slots  in  echelon 
[^12]  where  a'/X  “ 0.5  and  b'^X  ■ 0.05.  Our  computed  results 
are  compared  to  those  calculated  by  Borglottl  [1] . 


COUPLING  AMPLITUDE  I Su  I 


0.08 


0.06 


xPHRSr^°7MflILL0UX 
— COMPUTER  PROGRAM 


1^' 

L3 

a=a' 



5 3U  45  60  75 

ANGULAR  POSITION  0 (DEGREES) 


The  scattering  parameter  S^2  between  two  slots  where 
a/X  * a'/X  “ b/X  b'/X  ■ 0.6  and  S/X  « 0.9.  Our  computed 
results  are  compared  to  those  calculated  by  Mallloux  [3]. 


PHASE  OF  S«  (DEGREES) 
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Figure  8 shows  the  half-space  admittance  between  two  slots  In 
lis 

echelon  separation  distance  Is  varied  while  the  angular 

position  Is  held  fixed.  Three  different  slot  size  cases  are  shown. 

Figure  9 Is  the  same  as  Figure  8 except  that  the  angular  posi- 
tion of  one  slot  with  respect  to  the  other  Is  varied  while  the  separa- 
tion distance  Is  held  fixed. 

Figure  10  shows  the  coupled  power  (20  log  and  phase  of 

S^^(i  » 2,3,...,  49)  between  a driven  element  and  the  other  elements 
of  a 7 X 7 waveguide-fed  aperture  array  with  a rectangular  lattice. 

Figure  11  Is  the  same  as  Figure  10  except  that  the  lattice  Is 
now  Isosceles  triangular. 

VII.  DISCUSSION  AND  CONCLUSIONS 

A method  of  computing  the  mutual  coupling  for  an  array  of  uni- 
formly spaced  rectangular  waveguide-fed  apertures  radiating  Into  a half- 
space region  has  been  developed.  The  array  may  have  a rectangular  lattice 
or  an  Isosceles  triangular  lattice.  The  formulation  uses  one  expansion 
function  per  aperture  so  as  to  be  efficient  when  relatively  large  numbers 
of  apertures  are  concerned. 

For  calculating  the  half-space  self  admittance  for  a single 
aperture  an  eight-point  Gaussian  quadrature  numerical  Integration  [16] 

Is  performed  on  the  single  Integrals  In  (68) . The  same  method  of  solu- 
tion Is  used  for  calculating  the  half-space  mutual  admittance  for  aper- 
tures which  are  close  together  (centers  of  any  two  given  apertures  are 
separated  by  less  than  4a').  For  greater  aperture  separations  a six 


•-OT  X SOHW 
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hs 

Fig.  8.  The  half-space  admittance  between  two  slots  In  echelon 

where  A - (a'/X  - 1.0000,  b'/X  = 0.4761),  B - (a'/X  = 0.7500, 
b'/X  - 0.3571),  and  C = (a'/X  - 0.5000,  b'/X  = 0.2381). 


ANGULfiR  POSITION  0 (DEGREES) 


Fig.  9.  The  half-space  admittance  between  two  slots  in  echelon 
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point  double  numerical  Integration  (six  points  In  each  variable  for 
a total  of  36  points)  Is  performed  on  the  double  Integrals  In  (59). 

It  was  found  that  for  apertures  which  are  close  together,  the 
solution  Involving  the  single  Integrals  of  (68)  was  both  more  accurate 
and  efficient  than  the  solution  Involving  the  double  Integrals  of  (59) 

(an  eight  point  single  numerical  Integration  versus  a ten  point  double 
numerical  Integration).  As  the  separation  distance  between  apertures 
Increased,  the  complexity  of  the  solution  Involving  single  Integrals, 
specifically,  the  number  of  required  subdivisions,  was  not  necessary 
and  the  solution  Involving  double  Integrals  proved  Just  as  accurate 
and  more  efficient  (a  four  point  single  numerical  Integration  versus 
a six  point  double  numerical  Integration) . 

For  finding  the  unknown  coefficient  vector  V which  Is  required 
for  determining  the  scattering  matrix  [S],  the  computer  program  uses 
Zohar's  algorithm  [11]  for  the  linear  array  case  and  Slnnott's  algo- 
rithm [13]  for  the  rectangular  array  case.  Both  algorithms  utilize 

the  Toeplltz  symmetry  properties  of  the  admittance  matrix  [y'’®  + Y^®] . 

wc  lis 

Both  algorithms  are  more  efficient  than  Inversion  of  [Y  ^ -f  Y ] or  even 

a Gaussian  ellmlnatlon-LU  decomposition  algorithm  for  large  N.  (Zohar's 

2 

algorithm  requires  approximately  2N  multiplications  and  divisions, 

2 2 

Slnnott's  algorithm  requires  approximately  (n-fl)  multiplications 

and  divisions  where  n+1  Is  the  number  of  blocks  of  the  block-Toeplltz 

admittance  matrix  and  N Is  the  dimension  of  a block,  Gaussian  ellmlnatlon-LU 

P 

3 

decomposition  requires  N /3  multiplications  and  divisions,  and  a solution 

WC  Hs  3 

Involving  Inversion  of  [Y  ° + Y ] requires  N multiplications  and  divisions 
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% 

for  the  Inverse  and  another  multiplications  for  finding  V given  the 
Inverse.)  In  addition  to  the  efficiency  of  Zohar's  and  Slnnott's  algo-  | 

rlthms,  the  storage  requirements  are  considerably  less  than  Gaussian 
ellmlnatlon-LU  decomposition  or  Inversion.  (Zohar's  algorithm  requires 

W2  ll8 

the  storage  of  only  one  row  of  [Y  ® + Y ] while  Slnnott's  algorithm 
requires  storing  one  row  of  blocks  of  [y'^®  + Y^®].) 

If  very  large  arrays  are  encountered,  further  simplification 
of  the  program  could  be  made  by  using  the  present  formulation  for  aper- 
tures which  are  close  together,  and  using  a far-field  approximation  where 
the  field  Is  assumed  constant  over  the  aperture  for  the  rest  of  the 
apertures. 
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